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A  STUOY  OF  LONG  AEROSOL  INITIATED  LASER  INDUCED  AIR  BREAKDOWN  PLASMAS 


INTRODUCTION 

Partially  conducting,  reduced  density  channels  in  the  atmosphere  can  be 
produced  along  predictable  paths  by  high  energy  focussed  laser  beams.  These 
channels  are  interesting  for  at  least  two  reasons:  (1)  If  high  energy  particle 
beams  are  to  propagate  at  atmospheric  pressure,  they  will  have  to  propagate 
in,  and  they  will  produce,  reduced  density  channels,  (ii)  A  high  voltage 
electric  discharge  will  follow  the  path  of  a  reduced  density  channel.  Long, 
predictably  located,  electric  discharges,  guided  along  a  laser  produced 
channel,  provide  a  systematic  means  of  studying  lightning  phenomena.  The 
highly  conductive  afterglow  plasmas  of  these  precisely  located  electric 
discharges  have  been  used  as  radio  frequency  antennas  to  transmit  and  receive 
information.* 

For  several  years,  partially  conducting,  reduced  density  channels  in  the 
atmosphere  have  been  produced  at  NRL  using  both  Nd/glass  and  C02  lasers. 

These  lasers  produce  channels  that  are  reduced  in  density  to  about  40%  of 
ambient  density  and  have  conductivities  in  the  range  of  lxlO6  Hz.  An  electric 
discharge  guided  along  the  laser  produced  channel  further  reduces  its  density 
to  about  10%  of  ambient.  Longer  laser  produced  channels  are  possible  in 
aerosol  seeded  air  than  are  possible  in  clean  air.  Laser  produced  plasma 
channels  as  long  as  10  meters  have  been  created  In  air  seeded  with  aerosols 
produced  by  burning  black  powder. 

The  purpose  of  this  present  work  Is  to  Investigate  the  dependence  of  the 
length  of  the  laser  produced  plasma  channel  on  aerosol  density,  system  focal 
length,  and  laser  energy,  in  order  to  design  future  experiments  with  even 
longer  channels.  The  experimental  method  In  this  work  Is  to  correlate  the 
photographically  recorded  plasma  channel  lengths  with  the  above  mentioned 
parameters.  The  result  of  this  correlation  Is  a  determination  of  the  aerosol 
density  range  that  Is  useful  in  Increasing  the  plasma  channel  length  and  a 
determination  of  the  threshold  of  laser  Intensity  on  the  aerosol  particles 
that  produces  sufficient  ionization  to  cause  breakdown. 

Manuscript  approved  April  4,  1984. 


DESCRIPTION  OF  THE  EXPERIMENT 

The  apparatus  for  the  experiment  is  shown  schematically  in  Fig.  1.  An 
approximately  30  Joule,  60  nsec  FWHM,  1.06  micron,  laser  pulse,  after  passing 
through  power  and  energy  monitors,  is  turned  by  a  prism,  focussed  by  a  5  to  10 
meter  focal  length  lens,  and  finally  turned  again  by  a  prism  into  a  10  meter 
long  aerosol  chamber.  The  focus  of  the  optical  system  is  placed  at  the 
midpoint  of  the  chamber.  Open  shutter  cameras  photograph  the  laser  produced 
plasma-channel  to  record  its  length.  The  aerosol  density  is  known  as  a 
function  of  time  from  previous  measurements  which  used  an  optical  multichannel 
analyzer  to  measure  light  scattered  by  the  aerosols  from  an  argon  ion  laser 
beam. 

Following  in  this  section  is  a  description  of  the  experiment  chamber  and 
the  aerosol  generating  system,  a  description  of  the  measurement  of  the  plasma 
channel  length,  and  a  description  of  the  aerosol  density  measuring  system. 
Finally,  a  brief  description  Is  given  of  experiments  using  10  micron  diameter 
pollen  particles  instead  of  black  powder  smoke  as  the  breakdown  enhancing 
aerosol. 

A  schematic  diagram  of  the  experiment  chamber  is  shown  in  the  upper  part 
of  Fig.  1.  The  chamber  consists  of  six,  50  cm  diameter  by  3.05  meter  long 
sections  of  a  hard  paper  tube  used  in  the  building  industry  for  casting 
concrete  columns.  The  tubes  are  joined  Into  two  approximately  10  meter  long 
sections  located  side  by  side,  near  the  ceiling  of  a  large  room,  as  shown  in 
the  photograph  of  Fig.  2.  The  10  meter  sections  are  joined  together  at  each 
end  so  that  air  can  circulate  continuously  through  the  system.  The  upper  tube 
Is  the  laser-aerosol  interaction  region.  The  ends  of  this  tube  are  closed  by 
a  combination  turning  prism-window  on  one  end  and  a  Brewster  angle  blue  glass 
beam  dump  on  the  other.  A  ten  centimeter  high  acrylic  window  provides  visual 
access  along  almost  the  entire  length  of  the  upper  tube  for  photographing  the 
laser  produced  plasma  channel.  The  lower  tube  contains  a  fan  to  continuously 
circulate  the  aerosol -air  mixture,  and  a  port  connected  to  a  pump  that 
maintains  the  chamber  at  a  pressure  of  about  one  Torr  below  ambient 
pressure.  The  flow  rate  is  1.6  meters/second  or  one  revolution  per  12.4 
seconds.  One  end  of  the  lower  tube  is  sealed.  The  other  end  has  a  door  for 


introducing  the  black  powder  and,  just  inside  the  door,  an  asbestos  plate  on 
which  the  powder  is  burned.  To  fill  the  system,  1.5  grams  of  blackpowder  is 
ignited  with  a  hot  nichrome  wire. 

The  laser  induced  breakdown  experiments  are  undertaken  by  focussing  the 
pulsed  Nd/glass  laser  beam  into  the  center  of  the  aerosol  chamber.  Three 
parameters  are  varied:  laser  energy,  lens  f-number,  and  aerosol  density. 

Three  focussing  lenses  are  used,  5,  7.5,  and  10  meter  focal  lengths, 
corresponding  to  f-nunbers  of  138,  235,  and  360.  The  energy  of  the  laser 
pulses  is  varied  from  10  to  30  Joules.  The  energy  of  each  pulse  was  inferred 
from  its  peak  intensity,  as  measured  on  a  fast,  high  current,  ITT  F100 
photodiode.  Calorimeter  readings  and  laser  power  supply  settings  are  not  a 
good  indication  of  the  energy  in  the  laser  pulse  because  of  target 
Interactions  with  the  laser.  Pre-lasing  due  to  target  interaction.  In  this 
case  the  aerosols,  is  an  effect  that  causes  difficulty  only  with  the  short 
focal  length  experiments. 

The  raw  data  in  this  experiment  are  open  shutter  photographs  of  the  laser 
produced  plasma  channels  that  can  be  accurately  scaled  to  determine  their 
length.  The  cameras  are  located  7.65  meters  from  the  centerline  of  the 
experiment  chamber,  view  the  chamber  at  a  right  angle  to  Its  centerline,  and 
have  focal  lengths  of  50  and  162  mm.  To  calibrate  the  photographs,  a  scale 
with  10  cm  increments  Is  marked  on  the  chamber,  just  above  the  viewing 
window.  These  10  cm  marks  are  evident  in  Fig.  2.  The  plasma  channel  lengths 
can  easily  be  determined  within  2  cm  for  plasmas  produced  through  the  7.5  and 
10  meter  focal  length  lenses.  With  the  5  meter  focal  length  lens,  the  plasmas 
are  so  intense  that  light  scattered  from  the  windows  saturates  the  film, 
hiding  the  detail  of  the  spark.  Since  the  5  meter  focal  length  lens  produces 
only  relatively  short  plasmas,  there  Is  a  loss  of  accuracy  only  for  the  less 
interesting  data. 

The  aerosol  density  in  the  chamber  decreases  exponentially  from  the  time 
that  the  black  powder  is  ignited  until  It  reaches  the  density  of  aerosols  In 
laboratory  air.  From  experiment  to  experiment,  the  aerosol  density  Is  varied 
by  simply  waiting  until  the  desired  level  Is  reached.  The  Initial  aerosol 
density  is  determined  from  argon  Ion  laser  scattering  measurements  that  are 
described  in  the  next  section.  The  laser  attenuation  length  corresponding  to 
the  Initial  density  is  2.5  meters.  The  aerosol  density  decay  time  is  measured 
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to  be  about  8  minutes,  and  data  is  taken  from  10  seconds  after  black  powder 
ignition  to  10,000  seconds  after  ignition.  At  this  later  time,  the  aerosol 
density  *n  the  chamber  is  the  same  as  that  of  laboratory  air. 

Besi  black  powder  smoke  aerosols,  attempts  were  made  to  study  laser 
discharges  in  air  laden  with  ragweed  pollen.  Ragweed  pollen  grains  are 
approximately  spherical  and  about  10  microns  in  diameter:  over  ten  times  the 
diameter  of  black  powder  smoke  particles.  It  was  thought  that  the  larger 
particles  might  have  a  considerably  lower  threshold  of  laser  intensity  for 
producing  breakdown.  Unfortunately,  it  was  not  possible  to  properly  disperse 
the  pollen  and  the  agglomerates  settled  out  and  stuck  to  the  chamber  walls  too 
fast  to  make  a  density  determination.  The  few  discharges  that  were  made  with 
the  pollen  aerosol  were  not  as  long  as  corresponding  discharges  with  black 
powder  smoke. 

PROPERTIES  OF  THE  BLACK  POWDER  SMOKE  AEROSOL 

When  black  powder  (75%  KNO3,  15%  charcoal,  and  10%  sulphur)  is  burned,  a 
dense  white  smoke  is  evolved  that  consists  chiefly  of  potassium  sulphate, 
potassium  thiosulphate,  and  potassium  carbonate.  An  idea  of  the  complexity  of 
the  black  powder  reaction  is  given  by  Its  chemical  equation: 

74KN03  +  3 OS  +  16CgH20(charcoal)  56C02  +  14C0 
+  3CH4  +  2H2S  +  4H2  +  35N2  +  19K2C03  +  7K2S04 


+  2K2S  +  8K2S203  +  2KCNS*+  (NH^CC^ 

+  C  +  S  +  665  kcal/kg  .  (1) 

This  equation  Is  based  on  the  analysis  of  reaction  products  and  may  vary  under 
different  burning  conditions  and  compositions.  (See  Reference  2.)  Analysis  of 
this  reaction  shows  that  62%  of  the  initial  weight  of  the  black  powder  becomes 
particulate  material  after  burning. 

Microscopic  examination  of  smoke  particles  that  had  settled  on  a  glass 
surface  Indicates  a  particle  size  distribution  that  is  sharply  peaked  about 
one  micron  in  diameter.  However,  this  value  leads  to  conclusions  that  are 
either  unphysical  or  in  conflict  with  other  observations:  viz,  it  implies  that 


the  scattering  crossection  for  argon  ion  laser  light  is  ten  times  the  particle 
silhouette,  and  that  the  particles  should  have  about  equal  cross  sections  for 
scattering  both  .5145  micron  argon  ion  laser  light  and  1.06  micron  Nd/glass 
laser  light.  The  visual  measurement  of  the  particles,  therefore,  is 
disregarded  on  the  grounds  that  diffraction  effects  caused  errors  in  the 
measurements  or  that  the  particles  that  settle  out  are  not  representative  of 
the  total  aerosol  population.  Based  on  preliminary  electron  microscope 
measurements  and  on  the  observed  relative  magnitude  of  the  scattering  cross 
section  for  .5143  micron  and  1.06  micron  light  a  more  reasonable  value  for  the 
average  particle  diameter  is  0.5  micron,  and  this  is  the  assumed  value  for  the 
diameter  of  the  average  blackpowder  aerosol  particle. 

An  upper  limit  on  the  initial  aerosol  number  density  can  be  determined 
from  the  chamber  volume,  average  particle  diameter,  and  amount  of  powder 
burned.  For  a  chamber  volume  of  3.7  cubic  meters,  1.5  grams  of  powder,  and  an 
assumed  particle  diameter  of  0.5  microns  the  maximum  possible  initial  density 
is  1.5x10®  cm-3.  The  index  of  refraction  of  the  aerosol  particles  is  that  of 
the  constituent  potassium  salts,  about  1.5,  and  the  specific  gravity  of  the 
particles  is  about  2.5. 

DETERMINATION  OF  AEROSOL  DENSITY  BY  LIGHT  SCATTERING 

The  experiment  chamber  is  not  air  tight  and  is  maintained  slightly  below 
ambient  pressure.  A  result  of  this  constant  replacement  of  the  air  in  the 
chamber  is  that  Its  aerosol  density  decreases  exponentially  from  its  initial 
value  to  that  of  the  laboratory  air,  and  this  change  In  the  aerosol  content  is 
faster  than  any  changes  in  the  aerosol  size  distribution  caused  by  settling  of 
the  larger  particles.  To  determine  the  aerosol  density  in  the  chamber  as  a 
function  of  time,  light  scattered  out  of  the  beam  from  a  cw  argon  ion  laser  is 
recorded,  as  a  function  of  time,  from  the  time  of  ignition  of  the  black  powder 
to  the  time  that  the  aerosol  density  in  the  chamber  approaches  that  of  the 
ambient  air.  Figure  3  is  a  schematic  diagram  of  this  experimental 
arrangement. 

There  are  three  parameters  to  be  determined:  the  initial  aerosol  density, 
the  decay  rate  of  the  density,  and  the  total  scattering  cross  section  of  the 
aerosol  particles.  The  density  decay  rate  is  straight  forwardly  determined. 


as  it  is  the  decay  rate  for  the  scattered  light.  The  product  of  the  density 
and  the  total  cross  section  is  determined  at  a  single  instant  of  time  when  a 
max-mum  occurs  in  the  scattered  light  signal.  Since  the  cross  section  is 
assumed  to  be  a  constant  and  the  density  to  decrease  exponentially,  the 
density-cross  section  product  is  determined  for  all  time.  If  other 
measurements  can  determine  the  cross  section  or  the  initial  density,  the 
density  as  a  function  of  time  is  determined. 

Equation  (2)  is  the  expression  that  governs  the  size  of  the  light  signal 
detected  at  time,  t,  due  to  light  scattered  from  a  region  located  at  a 
position,  z,  in  a  volume  that  is  being  irradiated  by  a  beam  whose  intensity  at 
z=0  is  I0. 


I(z,t)  =  i0n(t)(do/dn)Ane”n(t)oz 

(2a) 

n(t)  *  n0e"Yt 

(2b) 

In  Equation  (2),  n(t)  is  the  exponentially  decaying  density  function,  y  is  its 
decay  rate,  and  (do/dn)AQ  is  the  scattering  cross  section  for  scattering  the 

laser  light  into  the  detector  optics.  The  exponential 

in  Eq.  (2a)  accounts 

for  the  attenuation  of  the  laser  beam  by  the  aerosols, 
maximum  at  the  value  of  t  given  by 

Equation  (2a)  has  a 

n(t)  a  z  =  1 

(3) 

Since  the  values  of  y  and  t,  at  the  maximum  of  I(z,t),  are  found  from  the 
scattered  light  record,  and  since  z  is  determined  by  the  experimental  set-up, 
nQo  is  determined  absolutely;  even  though  only  relative  intensity  measurements 
are  made. 

Figure  4  is  a  graph  of  the  scattering  signal  measurements  made  at  the 
three  axial  positions  along  the  aerosol  chamber  shown  in  Figure  3.  The 
measurements  are  made  at  three  minute  intervals  and  recorded  on  an  optical 
multichannel  analyzer.  The  data  points  are  compared  to  graphs  of  Eq.  (2a) 
with  y  equal  to  0.12  minute"*,  nQa  equal  to  0.4m"*,  and  z  equal  to  4.1,  5.25, 
and  6.3  meters.  The  value  of  n  a  is  determined  from  the  maxima  of  the  curves 


through  the  data  points  and  y  Is  determined  from  the  slope  of  the  curves  at 
late  time. 

If  the  initial  aerosol  density  from  the  previous  section  is  used  with  the 

measured  value  of  nrta  to  calculate  o,  the  result  Is  a  cross  section 
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of  2.6x10  cm  .  This  value  is  1.3  times  the  silhouette  of  a  0.5  micron 
diametwi'  particle,  a  rather  low  value  according  to  Mie  theory. ^  The  separate 
evaluation  of  nQ  and  o  is  interesting  because  it  is  a  check  on  the  efficiency 
of  the  aerosol  production.  On  the  other  hand,  it  is  only  academic,  since  the 
important  quantity  is  the  attenuation  of  the  laser  beam  by  the  aerosol,  and 
this  depends  on  the  product  of  nQ  and  o. 

RESULTS  OF  THE  EXPERIMENT 

The  organization  of  this  section  is  as  follows:  (i)  The  experimental  data 
is  presented.  (11)  The  more  meaningful  data  is  selected  for  further 
discussion.  (Ill)  The  useful  aerosol  density  range  is  determined,  (iv) 
Attenuation  of  the  laser  beam  by  the  aerosols  is  determined  to  be  negligible, 
(v)  From  the  known  intensity  distribution  of  a  focussed  Gaussian  beam,  the 
threshold  breakdown  intensity  Ibd,  is  determined. 

Appendix  I  Is  a  table  of  the  48  laser  shots  fired  into  the  aerosol 
chamber  for  this  experiment.  For  each  shot  the  following  are  given:  shot 
number,  optical  system  focal  length  and  f-number,  time  interval  between  firing 
the  laser  and  black  powder  Ignition,  laser  energy,  plasma  channel  length,  and 
the  channel  length  per  unit  laser  energy.  The  same  data  are  condensed  in 
Table  I.  This  table  is  a  matrix  Indexed  by  the  three  values  of  f-number  and 
four  time  intervals.  In  the  squares  are  the  minimum,  the  maximum,  and  the 
average  values  of  the  breakdown  channel  length  per  unit  energy  that  correspond 
to  laser  pulses  fired  with  the  specified  aperture  and  during  the  specified 
time  Interval. 

There  Is  a  correlation  In  the  data  of  Appendix  I  between  aperture  and 
laser  energy:  lower  energy  for  lower  f-number.  The  reason  for  this 
correlation  is  the  Interaction  between  the  laser  and  the  target  aerosols  that 
degrades  the  energy  In  the  Q-swItched  pulse  and  Is  more  severe  at  lower  f- 
n umbers. 


It  is  apparent  from  the  table  that  L'  (the  plasma  length  per  unit  laser 
energy)  is  greatest  for  the  aerosol  densities  that  occur  for  T  (the  time  after 
black  powder  ignition)  between  100  and  1400  seconds.  It  is  also  apparent  that 
the  values  of  L'  for  T  equal  to  10  seconds  are  inconsistent,  or  at  least  not 
simply  related  to  the  other  values  on  the  table.  The  reason  for  this 
inconsistency  is  that  the  aerosol  is  created  at  one  end  of  the  chamber  and 
requires  many  seconds  to  evenly  distribute  itself.  The  values  for  T  equal  to 
10  seconds  are  not  representative  of  a  uniform  aerosol  density,  as  are  the 
other  values.  The  values  of  L‘  for  T  greater  than  1400  seconds  are 
representative  of  a  laser  induced  breakdown  in  laboratory  air,  rather  than  an 
aerosol  enhanced  breakdown.  In  the  range  of  T  between  10  and  100  seconds,  the 
aerosol  density  is  high  enough  to  seriously  attenuate  the  laser  beam  before 
the  beam  reaches  focus.  It  is  the  region  of  T  that  is  between  100  and  1400 
seconds  and  its  corresponding  aerosol  density  range  that  produces  the  longest 
discharges.  This  is  the  range  that  is  interesting  and  it  is  the  data  in  this 
range  that  is  analyzed  further. 

The  first  thing  to  notice  about  the  discharge  lengths  in  this  100  to  1400 

second  region  is  that,  as  scattered  as  the  data  seems  to  be,  there  is  no 

correlation  between  the  time  of  the  discharge  and  its  length.  (See  Appendix 

I.)  Between  100  and  1400  seconds  after  the  black  powder  ignition  there  is  a 

factor  of  five  reduction  in  the  aerosol  density:  from  1.2x10®  to 
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2.4x10  cm  for  the  assumed  0.5  micron  particle  size.  If  the  Nd/glass  laser 
beam  were  attenuated  by  the  aerosol  as  the  argon  Ion  laser  is,  the  longer 
discharges  would  occur  at  later  times,  when  the  attenuation  is  less,  and  the 
values  of  Ibd  recorded  at  later  times  would  be  higher.  That  the  Nd/glass  • 
laser  beam  is  not  attenuated  is  consistent  with  the  assumed  0.5  micron 
particle  size  that  puts  the  scattering  process  in  the  regime  where  the 
scattering  cross  section  varies  Inversely  with  a  power  of  the  wavelength  that 
Is  between  two  and  four.  The  attenuation  length  for  1.06  micron  radiation  in 
this  density  range  varies,  at  worst,  from  13  to  64  meters,  much  longer  than 
the  five  meter  half  length  of  the  experiment  chamber,  and  the  corresponding 
variation  in  intensity  is,  at  most,  ±  15%. 

Figure  5  affirms  that  the  Nd/glass  laser  beam  is  not  strongly  attenuated 
by  the  aerosol.  The  figure  shows  four  photographs  of  laser  induced  breakdown 
plasmas  with  the  10  meter  focal  length  lens  system:  The  first  is  without  the 


aerosol  and  is  therefore  very  short.  It  does,  however,  locate  the  focal  point 
of  the  optical  system.  The  remaining  three  photographs  are  shots  with  the 
aerosol  present.  These  three  show  that  the  laser  induced  plasma  is  located 
symmetrically  about  the  focal  region.  If  the  laser  beam  were  attenuated,  the 
discharge  would  extend  further  toward  the  laser  from  focus  than  away  from  it. 


To  determine  the  threshold  intensity  for  laser  breakdown,  the  intensity 


distribution  of  the  focussed  Nd/glass  laser  beam  must  be  used.  The  TEM00  mode 


of  the  laser  produces  a  very  nearly  Gaussian  beam.  (Such  a  beam  has  an 
intensity  that  varies  radially  as  Iexp(-r2/r02) ,  where  rQ  is  the  nominal 
radius  of  the  beam  and  I  is  the  on-axis  intensity.  A  property  of  the  Gaussian 
function  is  that  Itt rQ2  is  the  total  power  of  the  beam.)  In  the  region  beyond 
a  focussing  lens,  the  on-axis  intensity  as  a  function  of  axial  position  is 
given  by  the  following  function: 


I(z)  -  I 


N2(l+N2)"1 


0  (z/RQ  -  N2(l+N2)"1)2  +  N2(l+N2)"2 


N  =  ira  /XR 


where  I0  is  the  intensity  at  the  lens,  RQ  is  the  radius  of  curvature  of  the 
beam  phase  front  immediately  after  the  lens,  N  Is  proportional  to  the  optical 
system  Fresnel  number,  X  is  the  laser  wavelength,  and  a  is  the  beam  radius  at 
the  lens.^  Table  II  contains  the  values  of  these  parameters  for  the  three 
optical  systems  that  were  used.  If  Equation  (4a)  is  evaluated  at  the  values 
of  z  and  I0  that  correspond  to  the  extremities  of  the  laser  induced  plasma  and 
the  energy  of  the  laser  pulse,  the  resulting  value  of  I ( z)  is  the  threshold 
Intensity  for  breakdown.  In  all  cases,  N  is  very  large  and  is  neglected. 

With  IQ  expressed  in  terms  of  the  laser  pulse  energy,  the  pulse  duration 
(60x10’^  seconds),  and  the  entrance  area,  and  with  R0/2a  expressed  as  the  f- 
number,  Eq.  (4a)  is  rearranged  as 


i  -  0.29  f/i  cibdy 


V 

y.  .  - 


S  V, 


where  1  is  the  total  length  (i.e.,  2z),  in  meters,  of  the  laser  induced 
plasma,  E  is  the  laser  pulse  energy  in  Joules,  and  Idd  is  the  breakdown 
intensity  threshold  in  GW/cm2.  The  breakdown  intensity  threshold  is  found 
from  the  data  by  (i)  plotting  the  observed  laser  induced  plasma  lengths 
against  their  corresponding  products  of  f -number  and  square  root  of  pulse 
energy,  (ii)  finding  the  straight  line  through  the  origin  that  best  fits  these 
points,  and  (iii)  evaluating  Eq.  (5)  along  this  line  and  solving  for  Ibd. 
Figure  (6)  shows  this  plot  and  the  "best  fit"  straight  line  to  the  f=235  and 
the  f*360  data.  The  slope  of  this  line  corresponds  to  a  laser  breakdown 
intensity  of  4.0x10®  W/cm2. 

In  summary,  the  density  of  black  powder  aerosols  that  extends  the  laser 
induced  plasma  lengths  Is  between  ,2xl05  and  2.4xl05  cm"3  (for  an  assumed 
particle  diameter  of  0.5  microns)  and  the  threshold  intensity  for  air 
breakdown  in  the  aerosol  atmosphere  is  4x10®  W/cm2. 

RESULTS  OF  A  VERY  HIGH  INTENSITY  LASER  EXPERIMENT 

Figure  (7)  is  a  photograph  of  a  10  meter  long  laser  induced  plasma  that 
was  made  with  a  200  Joule,  4  nanosecond  pulse  from  a  Nd/glass  laser.  The  beam 
Is  focussed  with  a  pair  of  lenses  that  produce  a  17  meter  focal  length  and  an 
f-number  of  213.  The  breakdown  intensity  threshold  is  estimated  to  be  4.5x10® 
W/cm2.  The  aerosol  is  produced  by  burning  9.4  grams  of  black  powder  in  a  long 
corridor.  The  resulting  maximum  aerosol  density  is  1.1x10"^  gram/cm3,  and  it 
produces  an  attenuation  length  of  11  meters  for  .6328  micron  He-Ne  laser 
radiation. 

Figure  (8)  is  a  schematic  diagram  of  the  experimental  layout.  Because 
the  experiment  is  conducted  in  a  long  corridor,  it  is  impossible  to  view  the 
entire  interaction  length  from  a  90  degree  position.  In  order  to  accurately 
measure  the  plasma  length,  lights,  at  one  meter  intervals  as  shown  in  Figures 
(7)  and  (8),  are  strung  along  the  left  wall  of  the  corridor.  The  position  of 
the  lights  relative  to  distance  along  the  optical  axis  is  calibrated  by  a 
photograph  showing  both  the  lights  and  a  string  with  one  meter  markers 
stretched  along  the  optical  axis.  In  Figure  (7),  the  intersection  between  the 
optical  axis  and  the  line  between  the  first  scale  mark  and  the  first  light  is 
nine  meters  from  the  last  focussing  lens.  Correspondingly,  the  intersection 


of  optical  axis  and  the  line  joining  the  thirteenth  light  and  scale  mark  is  21 
meters  from  the  lens. 


Two  things,  besides  its  remarkable  length,  differentiate  the  plasma  made 
by  the  200  Joule,  4  nanosecond  laser  pulse  from  the  plasma  made  by  the  30 
Joule,  60  nanosecond  laser  pulse:  (i)  The  former  is  almost  entirely  upstream 
from  the  laser  beam  focal  position,  while  the  latter  is  centered  about  the 
focal  position.  The  laser  induced  plasma  of  Figure  (7)  lies  eight  meters 
upstream  of  focus  and  two  meters  downstream,  (ii)  The  plasma  induced  by  the 
200  Joule,  4  nanosecond  laser  pulse  is  extremely  uniform  compared  to  the  lower 
power  laser  induced  plasma,  which  has  many  isolated  plasma  regions  along  its 
length.  The  reason  for  the  first  difference  is  probably  that  the  higher  power 
laser  ionizes  the  air  in  the  focal  region  to  such  a  high  conductivity  that  the 
beam  is  almost  entirely  attenuated.  The  reason  for  the  second  difference  may 
be  as  follows:  the  energy  absorbed  in  the  ionized  gas  around  a  single  aerosol 
particle  increases  with  the  third  or  fourth  power  of  time.5  In  the  long,  low 
power  laser  pulses,  each  individual  plasma  has  time  to  grow  and  to  absorb  a 
much  higher  energy  than  In  the  high  power  case.  These  large  plasma  volumes 
may  shadow  downstream  regions  of  the  beam,  causing  gaps  in  the  plasma 
brightness. 


DISCUSSION  AND  CALCULATION  OF  MAXIMUM  BREAKDOWN  LENGTHS 

A  short  focal  length  lens  converges  a  laser  beam  so  rapidly  that  the 
length  of  the  biconical  focal  region  that  has  an  Intensity  above  the  threshold 
for  producing  optical  breakdown  is  very  short.  As  the  focal  length  is 
increased,  the  breakdown  region  increases  in  length  until  a  point  is  reached 
where  diffraction  begins  to  reduce  beam  convergence.  Increasing  the  focal 
length  beyond  this  point  decreases  the  breakdown  region,  and,  if  the  focal 
length  is  increased  enough,  the  breakdown  region  will  vanish.  Somewhere 
between  these  extremes  is  an  optimal  optical  system,  which  is  the  subject  of 
this  section. 

The  breakdown  length  Is  a  function  of  several  parameters:  (i)  The  laser 
power,  P0,  measured  in  watts,  (ii)  The  threshold  for  optical  breakdown,  Ibd, 
which  was  determined  for  a  black  powder  aerosol  atmosphere  in  an  earlier 
section  to  be  4xl09  watts/cm2,  (ill)  The  radius  of  curvature  of  the  laser 
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beam  phase  surface  as  the  beam  leaves  the  last  optical  interface.  This 
length,  R0,  is  the  same  as  the  lens  focal  length  for  the  case  of  a  parallel 
beam  incident  on  a  simple  lens,  (iv)  The  wavelength  of  the  incident 
radiation,  X.  The  longer  the  wavelength,  the  more  important  are  the 
diffraction  effects,  (v)  The  Fresnel  number,  N,  defined  in  Equation  (4b), 
which  gives  a  measure  of  the  importance  of  diffraction.  The  smaller  the 
Fresnel  number,  the  more  important  are  diffraction  effects,  (vi)  The 
attenuation  constant  for  the  least  amount  of  aerosol  that  will  provide  the 
long  laser  induced  plasma.  This  number,  y,  was  determined  in  a  previous 
section  to  be  1.6xl0"4  cm"*.  Lastly,  for  practical  reasons,  only  the  shortest 
possible  systems  are  considered.  This  means  that  the  laser  intensity  on  the 
last  optical  component  must  be  as  high  as  practical,  i.e.,  equal  to  the  damage 
threshold,  Idam  =  0.2  x  109  watts/cm2. 

Figure  (9)  is  a  plot  of  the  intensity  of  a  focussed  Gaussian  beam  as  it 
leaves  the  lens,  converges  to  a  focal  region,  and  then  expands.  The 
horizontal  axis  is  distance  along  the  optical  axis,  normalized  to  a  quantity 
P0/Idamx»  that  is  the  1ens'  foca1  1en9th»  divided  by  N.  In  this  analysis,  the 
laser  power  and  wavelength  occur  only  as  normalization  factors  for  the  axial 
position,  so  that  all  actual  lengths  are  directly  proportional  to  laser  power 
and  inversely  proportional  to  laser  wavelength.  The  vertical  axis  is  the 
Gaussian  beam  intensity  normalized  to  the  optical  surface  damage  threshold 
intensity,  which  is  also  the  intensity  at  the  lens.  The  curves  are  plots  of 
Equation  (4a)  in  which  I0  is  set  equal  to  Idam,  and  wa2  in  Equation  (4b)  is 
set  equal  to  Po/I^am*  As  N  increases,  the  beam  behaves  increasingly  according 
to  geometrical  optics. 

The  breakdown  threshold  intensity  is  also  shown  in  Figure  (9).  Wherever 
the  intensity  curves  are  above  the  threshold  line,  laser  induced  breakdown 
occurs.  For  the  present  choice  of  a  breakdown  threshold  to  damage  threshold 
intensity  ratio  of  20,  breakdown  does  not  occur  for  N  less  than  5,  and  the 
maximum  breakdown  length  occurs  for  N  equal  to  six.  The  exact  positions  of 
the  end  points  of  the  laser  induced  plasma,  for  the  case  in  which  laser  beam 


attenuation  is  negligible,  is  obtained  by  setting  I(z')«Ibd  in  Equation  (4a) 
and  solving  for  z'.  The  result  is 


z'2>1  -  £  (B  t  /«B-8W) 

(6a) 

z'  -  zj  /ae-eW 

(6b) 

•  ■  Wbd  <  1.  s  ■  N2(1+N2)  _1 

(6c) 

where  Zi*  is  the  beginning  of  the  laser  induced  plasma,  z2'  is  the  end,  and 
z2'-Zi'  the  len9th  the  laser  induced  plasma.  Equation  (6b)  may  be 
solved  computationally  for  the  maximum  value  of  plasma  length.  For  a  damage 
threshold  of  0.2  GWcnT^  and  a  breakdown  threshold  of  4.0  GWcnf^,  the  maximum 
normalized  plasma  length  is  .05  for  a  value  of  N  equal  to  six.  The 
corresponding  value  of  the  maximum  laser  induced  plasma  length,  for  a  0.5GW, 
1.06  micron  laser  pulse,  is  .05  Pg/^am  x  wh1ch  equals  11.8  meters.  The  lens 
focal  length  for  this  case  is  39.3  meters. 

Figure  (9)  may  be  used  to  predict  laser  induced  plasma  lengths  up  to  12 
meters  since  the  attenuation  length  for  the  laser  beam  in  the  aerosol 
contaminated  air  is  at  least  64  meters.  No  assumed  values  of  any  experimental 
parameters  are  used  In  the  graph's  construction,  except  that  the  beam 
intensity  at  the  lens  is  kept  constant  as  N  varies.  This  graph  is  totally 
inadequate,  however,  for  predicting  the  100  times  longer  maximum  breakdown 
length  of  the  50GW  laser. 

Figure  (10)  is  a  set  of  Intensity  curves  that  may  be  used  to  determine 
laser  produced  plasma  lengths  In  cases  where  laser  beam  attenuation  cannot  be 
neglected.  This  figure  contains  an  assumed  value  of  o  of  .05  (Equation 
(6c)),  and  It  makes  use  of  the  experimental  procedure  of  allowing  the  beam  to 
converge  in  a  zero  attenuation  atmosphere  and  entering  the  attenuating  region 
at  a  value  of  intensity  equal  to  Ibd.  To  generate  the  curves  of  Figure  (10), 
each  curve  of  Figure  (9)  is  translated  to  the  left  by  an  amount,  z^,  (See 
Equation  (6a).)  such  that  z'*0  in  Figure  (10)  is  the  beginning  of  the  laser 
induced  plasma  for  each  value  of  N.  The  far  end  of  the  plasma  is  found  by 


equating  the  translated  intensity  function  including  attenuation  to  I6<j  and 
solving  for  a  value  of  z'  greater  than  zero.  The  equation  to  be  solved  is 


!dam  0  e 
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)2  +  B2/N2 


where  y  is  the  laser  beam  attenuation  constant  in  cm"*  and  the  axial  position 
normalization  factor,  6,  is  given  by 
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Figure  (10)  permits  Equation  (7a)  to  be  solved  graphically.  Equation  (7a)  is 
rearranged  as  shown  in  Equation  (7c): 


The  right  hand  side  of  this  equation  produces  the  intensity  curves  of  Figure 
(10).  The  left  hand  side  graphs  as  straight  lines  through  the  point  z'=u  and 
I(z' i/Idam*20;  For  a  9iven  attenuation  constant,  y»  the  straight 
line,  (■£  ey6z  )  is  drawn  on  the  graph  and  the  intersections  of  this  line  with 
the  intensity  curves  determine  the  normalized  position  of  the  far  end  of  the 
laser  induced  plasma  for  the  various  values  of  N.  The  actual  length  of  the 
breakdown  plasma  is  found  by  multiplying  the  normalized  length  by  6. 

The  exponential  function  corresponding  to  y  “(64  meters)'1  and  Po=0.5GW 
is  shown  in  Figure  (10).  It  intersects  the  N=6  Intensity  curve  at  z'*.0445. 
The  attenuation  of  the  laser  beam  by  the  aerosol  causes  a  10%  reduction  in  the 
length  of  an  approximately  11  meter  long  laser  induced  plasma.  The 
exponential  function  corresponding  to  the  the  same  attenuation  constant  and 
Po*50GW  is  also  shown.  The  breakdown  length  for  this  laser  is  maximized  by 
selecting  the  intensity  function  that  intersects  the  exponential  as  far  to  the 
right  as  possible.  The  Intensity  function  for  N*16  meets  this  requirement  and 
the  corresponding  values  of  z'  and  plasma  length  are  .015  and  354  meters.  The 
lens  for  this  case  has  a  focal  length  of  1470  meters. 
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The  values  for  laser  induced  plasma  lengths  computed  in  this  section 
represent  the  following  idealizations:  (i)  Lens  aberrations  are  neglected. 

(ii)  A  perfect  Gaussian  beam  is  assumed,  (iii)  Attenuation  by  "clean  air"  is 
neglected,  as  well  as  distortion  of  the  beam  by  atmospheric  fluctuations. 
Therefore  the  breakdown  plasma  lengths  computed  here  are  the  maximum  values 
possible  for  the  passive,  spherical  optical  systems  considered.  On  the  other 
hand,  these  computations  should  serve  as  starting  points  for  more 
sophisticated  computations  and  theories. 

CONCLUSIONS 

There  are  four  contributions  reported  here.  The  first  is  the  measurement 
of  the  laser  breakdown  threshold  for  1.06  micron  radiation  in  air  containing 
black  powder  smoke  aerosols.  The  value  of  the  threshold  is  measured  to  be 
4.0xl09  watts/cm2,  and  within  the  experimental  error  of  15%  is  the  same  for  60 
nanosecond  and  4  nanosecond  pulses.  The  second  contribution  is  an  estimate  of 
the  attenuation  constant  of  the  aerosol  for  the  1.06  micron  laser  beam  in  the 
regime  of  aerosol  densities  where  the  long  laser  produced  plasmas  are 
produced.  This  attenuation  constant  is  estimated  to  be  (64  meters)-1  at  the 
low  end  of  this  density  range.  The  third  contribution  is  a  method  based  on 
light  scattering  for  determining  the  attenuation  constant  of  an  aerosol  whose 
density  is  decreasing  exponentially  with  time.  This  method  has  the  advantage 
of  not  requiring  transmitted  light  measurements  or  absolute  intensity 
measurements.  The  attenuation  constant  is  determined  from  the  time  of 
occurrence  of  the  maximum  of  the  scattered  Intensity. 

The  fourth  contribution  is  a  prescription  for  maximizing  laser  induced 
plasma  lengths  under  idealized  conditions  and  a  graphical  system  for 
determining  such  plasma  lengths.  A  graphical  solution  finds  the  Fresnel 
number,  for  a  given  laser  breakdown  threshold  and  laser  damage  threshold  on 
the  focussing  optics,  that  maximizes  a  normalized  plasma  length.  The  actual 
laser  Induced  plasma  length  is  the  product  of  this  normalized  length,  the 
laser  power  divided  by  the  laser  wavelength,  and  the  damage  threshold.  The 
maximum  lengths  of  laser  induced  breakdown  plasmas  calculated  in  this  way  are 
of  the  order  of  hundreds  of  meters  for  lasers  existing  now  at  NRL.  These 
lengths  are  extremely  encouraging,  even  with  the  realization  that  the 


calculation  assumes  Ideal  conditions.  One  area  of  doubt  remains:  It  is  not 
known  for  sure  if  a  high  voltage  electric  discharge  will  follow  the  path  of  an 
optimized  laser  induced  plasma.  All  laser  guided  electric  discharge 
experiments  have  been  conducted  with  laser  pulses  that  deposit  an  energy 
density  much  higher  than  would  be  deposited  with  a  maximally  long  discharge. 
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TABLE  I 

This  table  Is  a  condensed  version  of  some  of  the  data  presented  In 
Appendix  I;  The  average  values  of  laser  Induced  plasma  length  per  unit  of 
laser  pulse  energy  are  Indexed  with  respect  to  time  Interval  and  f-number  to 
show  the  spread  In  the  data  the  maximum  and  minimum  values  are  also  shown. 


f-number 


T»10  sec 


10<T<100  sec 


100<T<1400  sec 


1400<T<104  cm 


APPENDIX  I 


This  appendix  contains  the  data  from  the  Nd/glass  laser  Induced  plasma 
experiment  that  Is  shown  In  Figures  (1)  and  (2).  The  column  titles  represent 
the  shot  number,  the  focal  length  of  the  laser  focussing  lens,  the  f-number 
for  that  lens,  the  time  Interval  after  the  Ignition  of  the  black  powder  that 
the  laser  was  fired,  the  energy  of  the  laser  pulse,  the  length  of  the  laser 


Induced  plasma,  and  the  le 
pulse. 

Shot  No.  F.L.  f-i 


of  the  plasma  per  unit 

Time  Energy 


'  sir 


energy  of  the  laser' 

Length  L1 

105  cm  3.54  cm/J 

305  9.93 

306  9.71 

264  13.54 

170  11.33 


Energy 


Length 


BEAM  DUMP  . 


BLACK  POWDER 
IGNITION 


■DIRECTION  OF  AIR  FLOW 


FOCUSSING  LENS 


5.9m 

70 


•EXHAUST 

FAN 


'  MANOMETER 


•  AIR  CIRCULATING  FAN 
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MULTICHANNEL  OPEN  SHUTTER 
ANALYZER  CAMERAS 


CALORIMETER 


LASER  BEAM 


FI  00 

PHOTODIODE 


BEAM 

SPUTTER 


Nd/GLASS 

LASER 


Figure  1. 


Schematic  diagram  of  the  laser  Induced  plasma  experiment.  The  Nd/glass 
laser  Is  monitored  and  focussed  Into  a  50  cm  diameter,  10  meter  long, 
controlled  atmosphere  tube.  The  length  of  the  resulting  laser  Induced 
plasma  Is  determined  from  open  shutter  camera  photographs.  The  black 
powder  aerosol  density  Is  Inferred  from  scattered  light  measured  by  an 
optical  multichannel  analyzer. 
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Figure  3.  Schematic  diagram  of  the  argon  ion  laser/aerosol  scattering  system  for 
measuring  the  product  of  aerosol  density  and  total  scattering  cross 
section.  The  laser  Is  incident  from  the  left  and  is  attenuated  by  the 
presence  of  the  aerosols,  as  shown  In  the  upper  part  of  the  diagram. 
Light  scattered  by  aerosols  in  small  volumes  around  Zj,  z2,  and  z^  are 
scattered  into  the  OMA.  The  scattered  light  Is  recorded  for  several 
density  decay  times. 
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CLEAN  AIR  BREAKDOWN 


AEROSOL  BREAKDOWN 


AEROSOL  BREAKDOWN 


AEROSOL  BREAKDOWN 


L  FOCUS  OF  LASER 
—  REGISTRATION  LINE 

Figure  5.  Photographs  of  laser  Induced  plasmas  In  the  experiment  chamber.  The  top 
plasma  Is  In  clean  air  and  Is  used  to  locate  the  laser  focal  point.  The 
three  other  plasmas  are  Initiated  on  aerosol  particles.  The  figure  shows 
that  the  laser  Induced  plasmas,  at  least  on  the  average,  extend  as  far 
beyond  focus  as  they  do  In  front  of  focus. 
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Figure  6. 


Plot  of  the  lengths  of  the  laser  Induced  plasmas  versus  the  product  of 
their  f-number  and  the  square  root  of  the  energy  of  their  laser  pulses. 
The  dashed  line  Is  the  best  straight  line  fit  to  the  data  points.  The 
slope  of  this  line  is  used  to  determine  the  Intensity  threshold  for  laser 
Induced  plasmas. 
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Figure  7.  Photograph  of  a  10  meter  long  laser  Induced  plasma  created  by  a  200 

Joule,  4  nanosecond,  laser  pulse.  Lights,  one  meter  apart,  are  located 
parallel  to  the  plasma  for  scaling  purposes. 


Figure  9.  Plot  of  the  intensity  along  the  optical  axis  of  a  focussed  Gaussian  beam 
for  various  values  of  Fresnel  number.  The  Intensity  is  normalized  to  the 
Intensity  at  the  lens,  and  the  axial  position  is  normalized  to  a  number 
that  is  related  to  focal  length.  The  breakdown  threshold  is  also  shown. 
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Figure  10.  Plots  of  the  intensity  along  the  optical  axis  of  a  focussed  Gaussian  beam 
for  various  values  of  Fresnel  numbers  and  with  the  position  coordinate  of 
each  curve  translated  such  that  I ( z'  =0) *Ib(1 .  The  laser  produced  plasma 
channel  for  all  curves  begins  at  z'*0.  The  straight  lines  represent  laser 
beam  attenuation,  and  the  intersections  of  the  attenuation  lines  and  the 
intensity  curves  define  the  far  ends  of  the  plasma  channels. 
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